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Abstract 

The structure of the inner edge of the accretion disk around a black hole can be altered, if the matter inside the 
marginally stable orbit is magnetically connected to the disk. In this case, a non-zero torque is exerted on its inner 
edge, and the accretion efficiency e can be much higher than that in the standard accretion disk model. We explore 
the radiation properties of an accretion disk at its sonic point around a black hole with a time-steady torque exerted 
on the inner edge of the disk. The local structure of the accretion flow at the sonic point is investigated in the frame 
of general relativity. It is found that the accretion flow will be optically thin at its sonic point for most cases, if 
the additional accretion efficiency 5e caused by the torque is as high as ~10 %. The results imply that the variable 
torque may trigger transitions of the flow between different accretion types. 
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1. Introduction 

Most work on accretion disks around black holes are on 
the assumption that there is no stress at the inner edge of the 
disk (Shakura & Sunyaev 1973). The inner edge of the disk 
should occur very close to the radius of the marginally stable 
orbit r„is, since the plunging matter in the region of unstable 
orbits rapidly becomes causally disconnected from the disk. 
Recently, Krolik (1999) questioned this assumption and sug¬ 
gested that the matter inside the marginally stable orbit can re¬ 
main magnetically connected to the disk. If this is the case, 
an additional torque would be exerted on the inner edge of the 
disk. Gammie (1999) has shown that the torque enhances the 
amount of energy released in the disk. The energy of the matter 
in the plunging region is extracted and then released in the disk. 
The extracted energy is transported outwards by the stresses in 
the disk. The torque alters the disk structure, and then the ra¬ 
diation spectrum of the disk (Agol & Krolik 2000, hereafter 
AKOO). In their calculations, AKOO assumed that the circular 
motion of the matter in the disk is Keplerian, and the released 
gravitational energy of the accreting matter is balanced by the 
radiation locally, i.e., no radial energy advection has been con¬ 
sidered in their calculations. 

In this work, we investigate the local structure of the accre¬ 
tion disk at the sonic point with a non-zero torque exerted on 
its inner edge. The radiation properties of the disk at its sonic 
point is explored. In our calculations, the radial energy advec¬ 
tion in the flow is considered, and the circular velocity of the 
flow is not limited to the Keplerian value. 


2. Model 

AKOO assumed that the circular motion of the matter in the 
disk is Keplerian instead of considering the radial force bal¬ 
ance, and the inner edge of the disk is therefore set at r„is 
manually. They further assumed that the generated energy is 
dissipated locally in the disk, i.e., the radial energy advection 
has not been considered in their work. The general relativistic 
effects are included in their calculations, as done by Novikov 
& Throne (1973). 

For a transonic accretion flow around a black hole, we take 
the radius of sonic point rg as the inner edge of the disk. We 
assume that a time-steady torque is exerted on the inner edge 
of the disk, as done by AKOO on the assumption of a Keplerian 
disk (and also see Paczynski (2000) and Afshordi & Paczynski 
(2002)). We include the angular and radial momentum equa¬ 
tions in our calculations and the circular velocity of the flow is 
not limited to the Keplerian value. Our calculations are in the 
frame of general relativity, and the radial energy advection in 
the flow is included in the energy equation. The local structure 
of the flow at the sonic point is therefore available. The torque 
may be variable and extended to a region of the disk outside 
the radius of sonic point. For simplicity, we restrict our calcu¬ 
lations to the simplest case as AKOO, i.e., a time-steady torque 
is exerted on the inner edge of the disk. 

The magnetic held plays an important role in angular mo¬ 
mentum transfer of the matter in the disk. The value of viscos¬ 
ity a can be determined by the magnetohydrodynamical pro¬ 
cesses taking place in the accretion flow (Hawley 2000). The 
value of a of the flow may vary with time and radius. The 
angular momentum transfer in the disk caused by any specific 
magnetohydrodynamical processes may be described approxi- 
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mately by the parameter a, though the precise value of a and 
its variation is unavailable unless the detailed magnetohydro- 
dynamical processes in the flow are included in the numerical 
simulations. The angular momentum transfer of the matter in 
the flow caused by the magnetohydrodynamical processes in 
the flow may correspond to a certain value of a, and the torque 
exerted on the inner edge of the disk is governed by the same 
magnetohydrodynamical processes. So, the values of a and 
the torque are determined simultaneously by the magnetohy¬ 
drodynamical processes taking place in the disk. In order to 
avoid being involved in the complicated magnetohydrodynam¬ 
ical processes taking place in the flow, we adopt a-viscosity 
and another independent parameter Siy to describe the angular 
momentum transfer of the matter in the flow and the torque 
on the inner edge of the disk, respectively. This may induce 
some inconsistency in our calculations, but we can explore the 
general radiation properties of the disk varying a for a given 
torque. This will not change our main conclusions on the ra¬ 
diation properties of the disk (see further discussion in Sect. 

5). 

For a disk with an additional torque exerted on its inner edge, 
the viscosity ly at the sonic point is 

( 1 ) 

where r'o(?'s) is the viscosity at the sonic point without an addi¬ 
tional torque, Siy is caused by the additional torque exerted on 
the inner edge of the disk. 

The main goal of this work is to explore the local structure 
and radiation properties of the disk at its sonic point. The radi¬ 
ation properties are mainly governed by the radiative processes 
and the energy dissipated in the flow. The a-description on the 
angular momentum transfer in the disk is therefore sufficiently 
good for present investigation, though the exact value of a is 
unknown. In this work, we focus on the local structure of the 
disk and the radiation properties at the sonic point of the flow. 
We assume that a-viscosity is valid in the flow at the sonic 
point. The variation of a along radius may not affect our main 
results of the radiation properties, since we focus on the local 
structure of the disk at the sonic point. 

3. Equations of accretion flows in Kerr geometry 

Many authors have included general relativistic effects in 
their models of accretion disks around black holes (e.g., 
Novikov & Throne 1973; Lu 1985; Abramowicz et al. 1996; 
Peitz & Appl 1997; Gammie & Popham 1998; Popham & 
Gammie 1998; Manmoto 2000; etc.). In this work, we mainly 
adopt the equations presented by Abramowicz et al. (1996, 
hereafter A96). One may refer to A96 for details. 

The metric of a Kerr black hole on the equatorial plane takes 
the form given by Novikov & Thorne (1973) {G = c = 1), 

ds'^ =—Lodt)'^+ ^dr'^ + dz'^, ( 2 ) 

A r‘‘ A 

where 

A = r‘^ -2Mr + a^, (3) 

A = r^ + r^a'^ + 2Mra^, (4) 


ijj = 


2Mar 

A 


(5) 


Here M is the black hole mass and a is the specific angular 
momentum of the Kerr black hole. 

The Lorentz factor 7 is defined by 


1 


( 6 ) 


where 


The quantities R and H are defined by 

42 

R^ = ——, 0 = 0 -a;. 

r^A 

The radial velocity component V is given by 

- - -= 7i;(’') = 


VT^ 

Combining Eqs. (6)-(9), we have 
, 1 -! 


72 = 


-C 2 


and 


V = 




\/l-02772 
The angular momentum conservation is 

M dl Id/ ,3/2A^/"7 ^dO\ ^ 

27rr dr r dr [ r^ dr ) ~ ’ 


2Trr 

where 

; = 7 


^3/2 \ 

r.3Al/2 ) 


(7) 

(8) 

(9) 

( 10 ) 


( 11 ) 


( 12 ) 


(13) 


is the specific angular momentum per unit mass, and the term 
of the angular momentum carried by the vertical flux of radia¬ 
tion has been neglected (A96). 

Integrating Eq. (12), we have 


M 

27r 


{l-k) 


-S^A3/2 


A3/273 dO 

r'^ dr ’ 


(14) 


where /q is the integral constant. We assume that an additional 
torque is exerted on the inner edge of the disk: r = rg, here rg is 
the radius of sonic point. At the sonic point, the viscosity i'{rs) 
is described by Eq. (1). The additional accretion efficiency 
caused by the torque is Se. We can obtain a relation at r = rg: 


SeM 

27rO 


-E(5z/A3/2 




(15) 


In the absence of an additional torque, the viscosity i'o{rs) at 
the sonic point is very low, namely zero-torque approximation. 
In this work, we consider the cases with significant additional 
torques on the inner edges of the disks, so Siy ^ vq is usually 
satisfied at the sonic point. For simplicity, we use an approxi¬ 
mation: ^{rs) ~ Siy in our calculations. We further assume that 
a-viscosity is valid even at the sonic point. The value of a at 
the sonic point may be different from that in the region outside 
the radius of sonic point. 


and 
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The radial component of the momentum conservation is 
V dV ^ 1 dP 

r T, dr ^ 


1 — V^ dr 
where 

yf = - 


(16) 


MA {n-n+)(n-n^) 




1 - 


(17) 


and P = 2Hp is the vertical integrated pressure. 

The angular velocities of the corotating and counterrotating 
Keplerian orbits are 

K-±^3/2 ±0^1/2- 

The surface energy generation rate and the cooling rate 
F~ are 

2 

(19) 


.^2 


= 7" -- 


dn 

dr 


and 


F- = 


l6aT^ 
3t ’ 


( 20 ) 


respectively. Here we only consider the optically thick case. 
The energy equation can be written as 

?adv _ _ jp— 


jpa 


The equation of state is 

1 ^4 , pkT 

P=i;aT H-, 

d pnin 


( 21 ) 


( 22 ) 


where the ratio of the specific heats of the gas 7 g = 5/3 has 
been adopted. The magnetic pressure has been neglected, since 
we restrict our calculations to the optically thick case. For the 
optically thick case, the radiation processes are not relevant to 
the magnetic pressure. 

Using Eq. (22), the entropy gradient can be calculated by 


TdS=- 


12 — —f} ] c?lnT — (4 — 3/9)c?lnp 


(23) 


P V 

where (3 = Pg/p. 

The advection cooling rate due to the radial motion of the 
gas is 


J.adv M ^dS 

27rr dr 


M p 
27rr p 




^-(4-30)^ 

dr 


dr 


The vertical force balance gives 


P 




(r^ + 0^)2 + 2Aa2 


pH^ 9 " ^3 (7-2 _|_ ^2^ _ ^q2 


(24) 


(25) 


Substituting Eqs. (22), (24) and (25) into Eq. (16), we can 
rewrite the radial motion equation as 


dr 


(26) 


where 


D = V- 


72-51/3-9/32 d 
56 - 45/3 - 3/32 U’ 


and 




-8 + 3/3 + 3/32 ^ 

'56-45/3- 3/32^"^ 


In / 7 ' 




+ 


16 - 12/3 2TTrF^^^ 


(27) 


^^ 72-51/3-9/32 ^ 

r 56 - 45/3 - 3/32 “ ’ 
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(r2+Q2)2 + 2Aa' 

(r2 + a?Y ~ 


56-45/3- 3/32 m ' 

The sonic point is defined by the condition 

D = A = 0, 


(28) 


(29) 


and we have the radial velocity of the flow at the sonic point; 


/ 72 - 51/3- 9/32 ^ 

V 56-45/3- 3/32 ) 


(30) 


The standard a-viscosity 

i0 = ^aCsH, (31) 

is adopted in this work. 

Substituting Eqs. (15), (19), (20) and (24) into Eqs. (28) and 
(29), the physical quantities ps and Tg of the flow at the sonic 
point are available, if the parameters M, M, a, Tg, a, Se and 
fls specified. 

The effective optical depth in the vertical direction is 


Teff — 0.5E-^ (/^es + 


(32) 


where Kes and ks are the electron scattering opacity and free- 
free opacity, respectively. 


4. Results 


In principle, the set of equations describing the accretion 
flow around a black hole can be integrated with suitable outer 
boundary conditions. The solution is required to satisfy the 
regularity condition (29) at the sonic point r = rg for a given 
Se. The specific angular momentum Is at the sonic point and the 
radius of sonic point Tg are simultaneously given by the numer¬ 
ical solution. In the case of Se = 0,ls — lo- The main numerical 
difficulties in solving the set of equations would be the partic¬ 
ular value of Is which should satisfy the regularity condition 
(29)(see A96 for the cases of Is ~ Iq). 

In this work, we search all possible solutions in the parame¬ 
ter plane Tg — Is instead of integrating the equations with outer 
boundary conditions numerically. Using this method, we can 
study the local properties of the disk at the sonic point. All 
possible solutions satisfying condition (29) can be found. The 
advantage of this approach is that all physical solutions satisfy¬ 
ing both the outer boundary conditions and the inner regularity 
condition (29) will be included in the solutions found in this 
way, though some solutions we And would be unphysical. 

As discussed in Sect. 3, we only consider the optically thick 
case (see Eq. (20)). So, all optically thick solutions should 
satisfy Ted > 1 at the sonic point. 

In this work, we focus on the geometrically thin accretion 
disk, all calculations are therefore performed for accretion rates 
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Fig. 1. The parameter plane rs — Is for all possible solutions with Fig. 3. The same as Fig. 1, but for a Schwarzschild black hole {a = 0). 

different values of (5e (labeled near the curves). The parameters: 

M = 10®Mq, a = 0.95 and M/M^dd — 0.001 are adopted. The 
parameters in the left side of the curves are for optically thick solu¬ 
tions, i.e., Teff > 1. The radius of the black hole horizon rh, radius of 
the minimal bound circular orbit Vmh and radius of the marginally sta¬ 
ble orbits Vms are marked in the figure. The upper panel is for ol = 1, 
while the lower panel is for cx. = 0.1. 



Fig. 4. The same as Fig. 2, but for a Schwarzschild black hole {a = 0). 


Fig. 2. The same as Fig. 1, but for M/MEdd = 0.1. 
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Fig. 5. The maximal radius of sonic point rs,inax as a function of 
6t for a = 0.95. The different values of the parameters are adopted: 
M/MEdd = 0.001, Q = 1 (solid line); M/Medd = 0.001, a = 0.1 
(dashed); M/MEdd = 0.1, a = 1 (dotted) and M/MEdd = 0.1, 
Q; = 0.1 (dot-dashed). 



Fig. 6. The same as Fig. 5, but for a Schwarzschild black hole (a = 0). 


M /MEdd <0.1. We found that all possible optically thick disk 
solutions obtained in this work do satisfy the thin disk approx¬ 
imation iT/r < 0.1 at the sonic point. The parameter planes 
Ts — Is for all possible optically thick solutions are present in 
Figs. 1-4, for different values of M/MEdd. Q;, and a. We found 
that the parameter plane Ts — h for all possible optically thick 
solutions is sensitive to these parameters. It is interesting to 
find that there is an upper limits on the radius of sonic point 
in the solution planes Tg — Is for given parameters. The upper 
limits on the radius of sonic point Vs^^ax for possible solutions 
are plotted in Figs. 5 and 6. 

5. Discussion of the results 

The physical, global solutions for optically thick flows are 
available by integrating a set of equations with suitable inner 
and outer boundary conditions. The problem is that the search 
for such a parameter plane for all physical global solutions 
with all possible different inner and outer boundary conditions 
would be very difficult. The parameter plane rg — 4 for phys¬ 
ical, global solutions will be included in the present parameter 
plane Tg — 4, but the global solution will set a more strict con¬ 
straint on the existence of optically thick flows. So, the conclu¬ 
sion is tight that no physical, optically thick disk solution at the 
sonic point will be present outside the parameter plane rg — 4 
obtained here, though the possibility cannot be ruled out that 
the flow is optically thin at the sonic point while it becomes 
optically thick at a large radius outside the sonic point. 

It is found that the location of the sonic point for an optically 
thick accretion flow moves towards the horizon of the hole with 
the increase of 5e. For high 5e, no optically thick accretion flow 
at the sonic point is present. A higher means that much en¬ 
ergy of the plunging matter inside the marginally stable orbit is 
extracted to the disk, and the disk is heated to a higher temper¬ 
ature. The free-free opacity Kg decreases with the increase of 
temperature. The disk then becomes optically thin. 

Comparing the results for different values of a, we And that 
the accretion flow with a low a will be optically thin even for a 
low i5e, i.e., the optically thick flows at the sonic point with low 
a are present only for low 5e cases. A high a makes the energy 
extracted from the plunging matter inside the inner edge of the 
disk transported outwards efficiently. Therefore, only a small 
fraction of the energy is released locally at the inner edge of 
the disk, and the flow can still remain optically thick for a rel¬ 
atively high 5e. For the lower a case, much energy is released 
locally and the disk is heated to a higher temperature, and the 
disk becomes optically thin at the sonic point. For a high ac¬ 
cretion rate m = M/MEdd, the surface density of the disk is 
so high that the disk is optically thick even for a relatively high 
5e (compare results in Figs 1-4 for different values of m). 

The sonic point of the accretion flow usually locates at the 
radius close to the marginally stable orbit, since the matter will 
plunge to the hole rapidly inside r^s- The radius of sonic 
point is in principle available by integrating the set of the equa¬ 
tions with suitable outer and inner boundary conditions. In the 
present work, we only consider the inner boundary condition at 
the sonic point, and only the range of the radius of sonic point 
is given. From Figs. 5 and 6, we And that 5e < 0.07 should be 
satisfied for an optically thick flow at the sonic point around a 
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Kerr black hole (a = 0.95), if rg > r^b is assumed (see Fig. 5). 
For a Schwarzschild black hole, Se < 0.02 should be satisfied 
(see Fig. 6). In our calculations, we found that the results are 
insensitive to the black hole mass. 

Recently, the XMM-Newton observation of the Seyfert 1 
galaxy MCG—6-30-15 reveals an extremely broad and red- 
shifted Fe Ka line (Wilms et al. 2001). The explanation on the 
observed spectrum requires the disk to have a very steep emis- 
sivity profile with an index around 4.3-5.0. A possible explana¬ 
tion is provided by the magnetic connection between the inner 
region of the disk and the plunging matter inside the inner edge 
of the disk (Wilms 2001; Krolik 1999; AKOO). Another pos¬ 
sible explanation is based on the magnetic connection between 
a rotating black hole and a disk (Wilms 2001; Li 2002a,b). 
However, it is also suggested that the origin of the broad Fe 
Ka line can be explained in the frame of an illuminated rela¬ 
tivistic accretion disk (Martocchia, Matt & Karas 2002). AKOO 
suggested that the locally generated surface energy flux scales 
as at large r in the limit of infinite efficiency or zero 

accretion rate in their model, rather than that scales as r~^ in 
the standard thin disk model. This is helpful to explain the 
required steep emissivity profile. In this work, we have inves¬ 
tigated the local structure and radiation properties of the disk 
at the sonic point. The global solution to this problem could 
be tested against the observed broad fluorescent Fe Ka line in 
MCG—6-30-15, which will be given in our future work. 

In this work, we find that the radiation properties of the disk 
at the sonic point are sensitive to the torque exerted on the in¬ 
ner edge of the disk. The torque may probably be variable (e.g., 
Armitage, Reynold & Chiang 2001; Hawley 2001; Hawley & 
Krolik 2001, 2002; Reynold & Armitage 2001; etc.). Our re¬ 
sults imply that the variable torque may trigger transitions of 
the flow between optically thick and optically thin accretion 
types. The global solutions are necessary to attack this prob¬ 
lem, which is beyond the scope of the present work, and it will 
be reported in our future work. 

In present calculations, we have assumed that the torque is 
exerted on the inner edge of the disk, as done by AKOO. If the 
torque is exerted on the extended region outside the inner edge 
of the disk, the disk can be optically thick at the sonic point 
even for a relatively higher 5e than that reported here. 
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